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Abstract 

A  non-linear  backward  equation  with  diffusive  terms  is  postulated  for 
the  probability  density  that  depends  on  the  Bohmian  quantum  potential. 

An  associated  nonlinear  Schrodinger  equation  is  also  introduced  and  ex¬ 
tension  of  the  analysis  to  several  particle  compounds  is  sketched  along 
with  the  implications  following  from  this  formalism  regarding  the  non¬ 
conservation  of  probability  in  the  non-equilibrium  regime.  Some  further 
conclusions  are  educed  with  respect  to  the  Onsager’s  relations  and  the 
generalized  optical  theorem. 

In  one  of  the  early  papers  on  wave  mechanics,  Darwin  mm  examined  an 
interesting  example  of  a  two-particle  quantum  interference  involving  a  set  up  of 
a  coherent  pair  of  particles.  A  few  years  later  interparticle  quantum  correlations 
between  two  bodies  became  notable  with  the  work  of  Einstein,  Podolsky  and 
Rosen[2l[3]  and  since  then  have  been  at  the  heart  of  quantum  mechanics.  The 
issue  of  non-separability  between  interacting  particles  had  been  particularly 
emphasized  by  Einstein  in  his  own  version  of  the  EPR  thought  experiment 
0.  Schrodinger  in  a  similar  vein  claimed  quite  succinctly  that  according  to 
quantum  theory  maximum  knowledge  of  a  composite  system  does  not  necessarily 
mean  fullest  knowledge  of  the  individual  particles  @133  ■  This  type  of  reasoning 
that  emphasizes  the  influence  of  the  particular  experimental  conditions  for  the 
available  type  of  possible  predictions  has  been  extensively  used  by  Bohr  in  his 
reply  to  Einstein-Podolsky-Rosen  [3]  and  is  best  illustrated  in  the  ingenious 
thought  experiments  devised  earlier  by  him  and  Rosenfeld  for  the  measurability 
of  fields  [S] .  More  recently  Mermin  in  particular  emphasized  the  physical  reality 
of  these  subsystem  correlations  in  comparison  to  the  correlata  underlying  them. 
From  a  rather  different  perspective,  the  quantum  potential  which  has  a  pre¬ 
eminent  role  in  Bohmian  dynamics,  is  associated  with  the  interrelatedness  in  the 
quantum  world  and  suggests  a  quality  of  irreducible  wholeness  in  the  domain 
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of  quantum  phenomena  mm-  Here  we  will  modify  the  continuity  equation 
for  the  probability  density  by  including  a  term  proportional  to  the  quantum 
potential  and  examine  some  of  its  consequences. 

Another  issue  that  needs  to  be  addressed  is  the  time  asymmetry  observed 
in  the  physical  world  in  contradistinction  to  the  time  symmetrical  form  of  the 
fundamental  equations  of  evolution  m-  The  resolution  of  the  contradiction  in 
classical  mechanics  is  usually  related  to  an  asymmetry  in  initial  conditions.  In 
quantum  mechanics  in  particular,  Schrodinger  equation  becomes  time  reversible 
if  we  make  the  identification  of  the  probability  amplitude  with  the  time  sym¬ 
metrical  complex  conjugate  one.  It  is  generally  appreciated  that  the  quantal 
irreversibility  is  the  result  of  initial  conditions  or  the  outcome  of  the  interaction 
of  a  few  particle  system  with  a  reservoir  but  recently  it  has  been  shown  experi¬ 
mentally  m  that  even  a  few  particles  may  exhibit  signs  of  irreversibility.  The 
collapse  of  the  wavefunction  on  the  other  hand  is  generally  accepted  to  impose 
a  genuinely  time  asymmetric  element  in  the  quantum  formalism. 

In  this  work  we  follow  a  different  approach,  by  demonstrating  that  the  quan¬ 
tum  potential  term  mentioned  above  may  be  employed  in  a  way  that  implies  that 
the  evolution  equation  for  the  probability  density  is  genuinely  time-asymmetric. 
It  deserves  notice  that  the  equation  that  will  be  proposed  depends  on  a  final 
value  condition  so  it  turns  out  to  be  backward  m  and  resembles  on  that  respect 
the  Kolmogorov  backward  equation  HE  the  adjoint  transport  equation  m  or 
the  Black-Scholes  one  known  from  mathematical  finance  mi.  Following  this  we 
extend  our  study  for  the  multiparticle  and  present  some  conclusions  related  to 
the  quantum  equilibrium  hypothesis,  Onsager’s  relations,  the  generalized  optical 
theorem  and  finally  we  recapitulate  our  findings  in  the  final  section. 

Generalized  continuity  equation  and  its  associated  non-linear  wave  equation. 
In  a  recent  paper  m  it  has  been  argued  that  quantum  theory  may  not  be  ex¬ 
tended  straightforwardly  to  the  macroscopic  domain.  This  could  lead  us  to  seek 
generalizations  of  the  present  quantum  formalism.  Since  standard  non-linear 
versions  of  the  Schrodinger  equation  ®  H3  seem  to  face  certain  difficulties 
mi  we  try  instead  to  generalize  the  Bohmian-de  Broglie  formulation.  Employ¬ 
ing  the  position  representation,  we  consider  initially  a  neutral,  spinless  particle 
of  mass  m,  assuming  that  the  time  varying  potential  is  V(t).  The  fundamen¬ 
tal  quantum  mechanical  equation  of  motion  for  the  single  particle  probability 
amplitude  reads  as 


(1) 


The  familiar  polar  substitution  T  =  RelS^h  gives  the  Hamilton- Jacobi-Madelung 
equation  when  we  equate  the  real  parts 
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We  may  also  define  the  quantum  potential  as 
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which  is  important  in  Bohm’s  version  of  the  de  Broglie-Bohm  mechanics,  al¬ 
though  it  is  possible  to  formulate  the  theory  without  any  reference  to  it.  Equat¬ 
ing  the  imaginary  parts  we  obtain 

flTD  1 

-  +  -(2VR,VS  +  RV2S)=0.  (4) 

Born’s  rule  which  at  least  in  text-book  quantum  mechanics  is  considered  a  sep¬ 
arate  axiom  of  quantum  physics  provides  a  conceptual  and  calculational  bridge 
between  theory  and  experimental  results.  Taking  into  account  then  the  identi¬ 
fication  p  =  u 2  for  the  probability  density  and  multiplying  both  sides  with  2 u 
yields  the  familiar  continuity  equation 


dp 
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From  (2)  and  (4)  we  may  determine  in  principle  the  action  S  and  the  amplitude 
R.  We  seek  to  modify  the  previous  expression  (4)  including  a  diffusion  term  and 
a  linear  term  proportional  to  the  quantum  potential  as  follows 


ft +  b  <2V“ ' vs + “v3s>  =  {~b)  (y3"  -  tt“)  ■  <6> 

This  formula  may  be  considered  as  a  non-equilibrium  extension  of  (4)  in  the 
sense  that  as  u  — >  R  the  former  reduces  to  the  familiar  equation  that  holds  at 
quantum  equilibrium.  The  given  expression  remains  time-asymmetrical  which 
immediately  implies  non-unitary  evolution  (the  inverse  proposition  is  not  true 
nevertheless).  The  negative  coefficient  suggests  a  backward  evolution  for  u(t). 
We  conclude  then  that  not  only  the  evolution  equation  is  not  time-invariant 
but  that  also  there  is  time  asymmetry  in  the  imposition  of  terminal  instead  of 
initial  conditions.  An  element  then  of  finality  is  introduced  in  (6)  since  it  relies 
on  terminal  values  evolving  backwards  in  time.  The  homogeneity  condition 
is  satisfied  since  for  an  arbitrary  solution  u  and  n  a  real  number  it  follows 
immediately  that  ku  is  also  a  solution.  Hence  the  rescaling  transformation  u'  — > 
ku  constitutes  a  similarity  transformation  that  leaves  the  equation  invariant.  An 
imaginary  quantum  potential  would  give  rise  to  dissipative  terms  that  do  not 
preserve  unitarity  in  the  general  case.  At  the  same  time  we  must  ensure  that 
the  expression  reduces  to  the  standard  form  at  equilibrium.  This  suggests  an 
extension  of  (1)  written  as 
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where  $  =  uelS^h  and  the  wavefunction  T  =  RelS^h  is  governed  by  (1)  with 
appropriate  final  conditions.  The  existence  of  a  distinct  wave  function  with 
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the  same  phase  as  HI  but  different  real  amplitude  is  reminiscent  of  the  double 
solution  theory  put  forward  by  de  Broglie  [22],  but  it  must  be  noted  that  in  (7) 
$  propagates  backward  in  time  in  the  configuration  space  in  the  same  manner 
as  'k,  so  it  is  not  a  singularity  in  the  real  three  dimensional  space  as  de  Broglie 
envisaged.  Dynamical  terms  of  the  form  given  in  (7)  have  also  been  used  in 
the  past  in  the  real  but  not  in  the  imaginary  part  of  the  Hamiltonian  123].  It 
must  be  underscored  that  when  the  terminal  condition  is  fulfilled,  which  means 
u(t  — >  oo )  =  R{t  — >  oo )  so  <f>(f  — >  oo )  =  H >(t  — ►  oo )  the  superposition  principle 
holds  and  we  recover  (1)  from  (7)  even  though  the  coefficient  h/m  is  obviously 
different  from  zero.  Hence  the  magnitude  of  the  non-linearity  coefficient  may 
be  non-zero  nevertheless  there  is  an  appropriate  limit  that  we  recover  (1).  On 
the  other  hand  a  nonlinear  Schrodinger  as  in  [19]  reduces  to  the  linear  one  only 
when  the  nonlinearity  constant  becomes  zero.  Finally,  (7)  is  time  irreversible 
even  when  — >  $*  and  respects  the  homogeneity  condition  |20|  which  says  that 
if  $  is  a  solution  of  (7)  so  is  — >  Z<h,  where  Z  a  constant  complex  number. 
Based  on  (7)  we  can  present  extensions  of  standar  quantum  formulas.  If  we 
write  for  example  the  expectation  value  of  an  operator  A  as 

(i)  =  J  <f>*A$d3x,  (8) 

its  time  development  is  given  by 
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where  H  the  Hamiltonian  and  of  course  the  first  two  terms  make  up  the  familiar 
expression  in  the  Schrodinger  picture  whereas  the  last  term  reflects  dissipative 
effects.  According  to  (6)  and  taking  into  account  that  p  =  |4>|2,  the  correspond¬ 
ing  expression  for  the  probability  density  is 
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(10) 


The  left  part  is  familiar  from  the  continuity  equation  but  the  right  non  zero  part 
implies  non-unitary  evolution.  The  existence  of  these  terms  suggests  a  possible 
interplay  between  spatial-temporal  coherence  and  diffusion  phenomena.  Diffu¬ 
sion  terms  in  the  continuity  equation  have  been  also  been  added  in  [241 125]  (see 
also  PB]  [27]).  In  addition  equation  (10)  remains  invariant  when  R  — >  kR  and 
p  — >  A p.  It  must  be  noted  though  that  the  diffusion  constant  in  our  proposal 
is  negative,  namely  —h/m  as  we  mentioned.  Without  the  the  nonlinear  and 
the  quantum  potential  term,  (10)  is  a  non-equilibrium  backward  Fokker-Planck 
type  equation  that  conserves  probability  locally.  The  inclusion  of  the  quantum 
potential  and  nonlinear  terms  destroy  this  property  and  is  responsible  for  its 
non-conservation.  It  is  evident  that  superluminal  communication  is  precluded 
by  the  very  fact  that  the  relevant  equation  of  motion  for  the  probability  density 
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depends  on  the  imposition  of  a  final  condition.  This  preferential  status  of  ter¬ 
minal  conditions  instead  of  initial  means  that  we  could  use  (6), (7),  or  (10)  to 
make  retrodictions  but  not  predictions.  It  may  be  commented  that  this  special 
importance  of  final  conditions  can  be  clearly  illustrated  with  the  help  of  the  well- 
known  thought  delayed  choice  experiment  128]  and  even  the  Bohrian  emission  of 
photons  during  atomic  transitions  in  the  old  quantum  theory  era.  At  the  equi¬ 
librium  limit  these  reduce  to  the  more  familiar  Bohmian  equations  so  they  may 
been  seen  as  generalizations  of  the  latter  in  the  non-equilibrium  domain.  The 
presence  of  the  diffusion  terms  would  tend  to  make  u  smoother  but  the  quantum 
potential  driving  term  and  especially  the  highly  nonlinear  midterm  could  have 
the  opposite  effect.  On  physical  grounds  we  would  expect  irregularities  in  p(x,  t) 
to  become  less  smooth  with  backward  time  evolution  which  may  cause  problems 
for  the  behaviour  of  the  solutions  like  arbitrary  rapid  growth  in  finite  time  [14] . 
The  presence  then  of  the  nonlinear  terms  is  of  special  importance,  since  if  they 
tend  to  create  irregularities  as  we  propagate  backwards  in  time  and  dominate 
over  the  decay  term  that  will  decide  in  favour  of  the  ill-posedeness  of  the  prob¬ 
lem  and  would  demand  sophisticated  regularization  techniques.  Obviously,  we 
could  also  include  advanced  vector  potentials  and  even  noise  terms  in  (10)  [29]. 
The  latter  ones  require  special  care  though  because  of  certain  mathematical 
subtleties  involved  mi- 

it  must  be  emphasized  though  that  there  are  no  non-intersecting  particle 
trajectories  apart  from  the  case  of  equilibrium  when  (8)  becomes  the  familiar 
density  current  equation  and  the  coefficient  sets  time  scale  equal  to  ft/m.  We 
now  seek  a  guidance  equation  modifying  the  de  Broglie-Bohm  one.  Motivated 
by  (11)  we  write  a  formula  for  the  velocity  v  of  a  single  particle  that  would 
accompany  (7)  as 


(V|$|2)2  2]$|  2V2yW'| 

2|<h|2  +  ^  )  • 

(11) 

The  first  term  in  the  above  formula  is  the  velocity  term  — ,  the  second  one 
proportional  to  VZn|<f>|2  is  an  osmotic  type  contribution  to  velocity  according  to 
[TUI  jSTJ  02]  (note  the  negative  sign  of  the  diffusion  coefficient)  and  the  third  inte¬ 
gral  term  expresses  the  non-local  character  of  this  generalized  guidance  equation 
and  clearly  violates  the  uniqueness  propery  of  trajectories.  Since  the  last  two 
term  cancel  each  other  when  <!>  =  we  recover  the  standard  guidance  equa¬ 
tion.  Having  established  the  basic  elements  of  the  formalism  (relevant  evolution 
equations)  we  proceed  by  examining  some  further  developments  and  examples. 

It  is  interesting  to  note,  as  it  is  clearly  seen  in  (10)  or  (11),  that  even  when 
the  gradient  of  the  action  is  zero,  the  probability  density  may  be  nonstationary 
which  implies  the  possibility  of  backflow  effects.  The  inverse  is  also  true,  (10) 
allows  a  steady  state  density  even  when  the  the  momentum  is  non-zero.  Since 
ft/m  is  very  small  we  expect  that  these  deviations  from  (5)  to  be  very  weak. 
Seeking  a  similarity  solution  we  proceed  by  defining  the  scaling  anzatz  r/  = 
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and  obtain  a  nonlinear  ordinary  differential  equation  which  still  may  be  not 
amenable  to  analytic  solution.  Instead  we  could  make  the  substitution  p  =  e2h 
to  (10)  when  the  spatial  derivative  of  the  phase  is  considered  negligible  and  find 


dh 

~dt 


V2h  +  (Vh)2  +  -TQ 
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(12) 


The  above  resembles  the  celebrated  Kardar-Parisi-Zhang  equation  [53]  that  has 
many  applications  in  the  field  of  surface  growth.  The  crucial  difference  is  that 
it  propagates  backward  in  time  and  the  last  quantum  potential  driving  term 
should  be  considered  deterministic  in  general  instead  of  stochastic. 

Non-conservation  of  total  probability  and  entropy  production  rate.  The  rate 
of  change  of  the  total  probability  is  expressed  as 


dP  d  J  d3xp  [  3  dp  (  h 

dt  dt  J  X  dt  \  m 


2  '  4  p 


V2R 

—R~f 


d3x. 


.  {13) 

The  second  term  makes  apparent  a  connection  to  Fisher  information,  which  is 
a  measure  of  system  disorder  [34j  [35] .  As  p  aprroaches  R  the  rate  of  change  of 
the  total  probability  becomes  zero.  Making  use  of  the  identities  =  NlnR 
and  =  V2lnR  +  (VlnR)2  and  allowing  p  — >  R2  +  Sp,  the  variation  of  the 
total  probability  rate  51  may  be  written  as 


"  =  (“£)  /  ^ 

Differentiating  over  time  the  Gibbs  entropy  G  =  —  J  d3xp\np  which  relies  on 
the  spatial  distribution  of  the  probability  density  gives  the  following  rate  for  the 
production  of  entropy 


Qv2  -  VlnRV  -  V2lnR^J  Sp.  (14) 


2  \/2R  1 


,(15) 

The  second  integral  remains  unaffected  when  time  is  reversed  and  can  be  iden¬ 
tified  as  a  measure  of  entropy  production  E  and  a  sign  of  system  irreversibility. 
The  other  term  in  (15)  is  associated  with  entropy  flow  that  follows  from  (5). 
Again  near  the  linear  approximation  regime  p  — >  R2+5p  the  variation  of  entropy 
production  is  given  by 


dG  f  ,3  /  VS\  h  f  3  ,  J 

— —  =  —  d3x  Vp -  1-1 - /  d3x(  1  +  Inp) 

dt  J  \  m  J  m  J 


2  F  4  p 


5Y>  =  VFiUjVFj .  (16) 

ij 

In  the  last  equality,  are  generalized  thermodynamical  forces  near  equilibrium 
|36j  and  Lij  the  corresponding  Onsager  coefficients.  Since  microscopic  time 
reversibility  and  the  resulting  detailed  balance  do  not  hold  it  follows  according 
to  (15), (16)  that  these  coefficients  may  not  be  symmetrical. 
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Quantum  scattering  and  the  generalized  optical  theorem.  Up  to  this  point  we 
have  only  considered  time  dependent  evolution  equations.  It  is  possible  never¬ 
theless  to  make  use  of  a  time-independent  version  of  (7)  invoking  its  imaginary 
potential.  The  time-independent  variant  of  this  nonlinear  wave  equation  for  a 
particle  with  energy  E,  mass  m  and  wave- vector  k  is  written  as 


(V2  +  k2  -  Ur) 


(vWWW 

v  vww 


V\*W  ) 


$(r)  =  0 


(17) 


where  k  =,  U( r)  =  ^™V(r)  and  V(r)  a  central  potential.  When  <f>(r)  =  ’F(r) 
the  formal  solution  of  (17)  can  be  written  as 


\&(r)  =  4f0(r)  +  J  Gk(r,r')U(r')^/(r')d3x  (18) 

where  Gk(r,r')  the  corresponding  Green  function  and  'l'o(r)  a  solution  of  the 
homogeneous  equation  (V2  +  fc2)\Ho(r)  =0  [37].  Due  to  the  nonlinear  character 
of  (17)  it  is  not  possible  to  write  a  relation  similar  to  (18)  for  $(r).  We  can 
proceed  though  by  assuming  that  at  large  distances  from  the  potential  region 
the  solution  $(r)  has  the  asymptotic  form 

‘hk  — ►  eikr  +  ——f(k,k').  (19) 

r 

Following  the  analysis  by  Glauber  [38]  [33  for  two  scattering  vectors  k,  k'  (17) 
produces  the  identities 


$_k'V2$k-$k'V2^k  =  0  (20) 

and 


V  •  [$i,V$k  -  *kV$£,]  +  2 % 


vFF  ) 
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(21) 


Integrating  (20)  and  using  Green’s  theorem  over  the  volume  of  a  very  large 
sphere  gives  /( k',  k)  =  /(— k',  — k)  which  expresses  the  reversibility  of  the  scat¬ 
tering  processes/amplitudes  for  two  directions  [38] .  Furthermore,  if  the  scatter¬ 
ing  potential  has  inversion  symmetry  then  the  scattering  amplitude  is  symmetric 
/( k',  k)  =  /( k,  k')  which  is  deary  related  to  the  principle  of  detailed  balance 
[38].  This  conclusion  nevertheless  does  not  contradict  our  previous  discussion 
which  relied  on  a  time-dependent  evolution  equation.  Following  the  same  pro¬ 
cedure  for  (21)  we  find 


i[/(k',k)  -  /*  (k,k')] = Ay  /*(k",k)/(k",k)df!k'' 
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Setting  k  =  k'  gives 


Im/(k,k') 
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(23) 


The  integrand  in  the  right  hand  of  (22)  is  easily  identified  with  the  divergence  of 
the  current  density  j  =  |<I>k|2^f  ■  Writing  <f>k  =  esWvj/k  (since  they  both  have 
the  same  phase)  where  g( r)  an  unknown  real  function  to  be  determined  yields 
the  integral  equation 


47rlm/(k,  k')  =  fccr 


V2.g(r)  +  Vfl(r)  •  +  (V.g(r))2 


e2s(r)  |\J/k|2<i3cc. 


(24) 

It  is  clear  that  we  recover  the  optical  theorem  for  elastic  scattering  when  $k  = 
Tk  which  it  is  known  that  expresses  the  particle  conservation  requirement  but 
more  generally  deviations  may  be  observed  according  to  (23)  or  (24).  When 
the  above  quantity  is  positive  attenuation  of  the  incident  beam  is  implied.  It 
is  interesting  to  note  that  the  presence  of  the  imaginary  potential  in  (17)  intro¬ 
duces  non-unitarity  since  the  total  particle  number  is  not  conserved,  but  does 
not  affect  the  reversibility  of  the  scattering  amplitude  [38]. 

Many  particle  case  and  conclusions.  It  is  straightforward  to  extend  our 
previous  results  for  the  case  of  a  compound  of  several  interacting  particles  of 
equal  masses  m 


The  presence  of  the  many  particle  quantum  potential  in  the  continuity  equa¬ 
tion  signifies  an  important  feature,  namely  the  dependence  of  the  multiparticle 
density  on  a  quantity  that  exhibits  non-local  characteristics  (but  since  (25)  is 
backward  in  time  there  is  no  possibility  of  non-local  signalling).  The  intercon¬ 
nectedness  of  the  system  constituents  is  expressed  then  through  the  form  of  this 
non-local  quantum  potential.  The  constituent  particles  then  are  reidentified 
and  incorporated  into  a  substantial  whole  unified  by  a  single  unifying  principle 
and  cannot  been  considered  as  distinct  aggregated  components.  This  means 
that  the  component  parts  even  though  are  retrievable  in  principle,  they  are  not 
intrinsically  unaffected  by  their  union  13300  HD-  Similar  considerations  in  a 
broader  sense  have  also  been  expressed  various  authors.  Whitehead  for  example 
has  spoken  about  the  concept  of  a  dynamic  and  contextual  organic  whole  [T21. 
as  of  course  Bohm  and  Hiley  [43]  01]  already  mentioned. 

In  summary,  in  the  present  genuinely  multiparticle  formulation  ontological 
priority  is  given  not  to  the  monadic  particle  but  to  the  emergent  several-particle 
composite.  We  conclude  that  the  presence  of  the  many  particle  quantum  po¬ 
tential  in  the  modified  backward  continuity  equation  suggests  non-separability 
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and  illustrates  the  irreducible  nature  of  a  multi-particle  compound.  It  is  rea¬ 
sonable  to  deduce  then  that  we  cannot  examine  each  particle  singulatim  neither 
we  should  consider  the  quantal  multiparticle  as  a  mere  mereological  aggregate 
of  distinct  co-present  entities  accidentally  unified,  but  rather  as  a  substantial 
unified  whole  that  constitutes  the  production  of  a  new  physical  reality.  This 
particular  formalism  also  ensures  that  the  time  anisotropy  inferred  from  the 
evolution  equation  and  its  dependence  on  terminal  values  is  nomological  and 
intrinsic  in  nature  rather  than  contingent  and  accidental. 
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